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Abstract 

Analytical expressions for the non-relativistic and relativistic Sunyaev- 



. Zel'dovich effect (SZE) are derived by means of suitable convolution in- 

I 1 tegrals. The establishment of these expressions is based on the fact that 

the SZE disturbed spectrum, at high frequencies, possesses the form of a 
Laplace transform of the single line distortion profile (structure factor). 
^ ' Implications of this description of the SZE related to light scattering in 

optically thin plasmas are discussed. 

^ ! 1 Introduction 



Distortions to the cosmic microwave background radiation (CMBR) spectrum 
arise from the interaction between the radiation photons and electrons present 
in large structures such as the hot intracluster gas now known to exist in the 
universe. Such distortions are only a very small effect changing the brightness 
of the spectrum by a figure of the order of 0.1 percent. This effect, excluding 
the proper motion of the cluster, is now called the thermal Sunyaev-Zel'dovich 
effect and its detection is at present a relatively feasible task due to the 
modern observational techniques available. Its main interest lies on the fact 
that it provides information to determine important cosmological parameters 
such as Hubble's constant and the baryonic density H]-®- 

The kinetic equation first used to describe the SZE was one derived by 
Kompaneets back in 1956 [3] 0-[E| This approach implies a photon diffusion 
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description of the effect that works basically when the electrons present in the 
hot intracluster gas are non-relativistic. Many authors, including Sunyaev and 
Zel'dovich themselves, were very reluctant in accepting a diffusive mechanism 
as the underlying phenomena responsible for the spectrum distortion. Later 
on, Rephaeli and others [5] computed the distorted spectrum by considering 
Compton scattering off relativistic electrons. Those works show that, at high 
electron temperatures, the distortion curves are significantly modified. Useful 
analytic expressions that describe the relativistic SZE can be easily found in the 
literature [TO]. 

One possible physical picture of the effect is that of an absortion-emission 
process in which a few photons happen to be captured by electrons in the opti- 
cally thin gas. An electron moving with a given thermal velocity emits (scatters) 
a photon with a certain incoming frequency v Q and outgoing frequency v. The 
line breath of this process is readily calculated from kinetic theory taking into 
account that the media in which the process takes place has a small optical 
depth directly related to the Compton parameter y. When the resulting ex- 
pression, which will be called structure factor, is convoluted with the incoming 
flux of photons obtained from Planck's distribution, one easily obtains the dis- 
turbed spectra. For the non-relativistic SZE, a Gaussian structure factor has 
been successfully established JT], but the obtention of a simple analytic rela- 
tivistic structure factor is a more complicated task J2]- Nevertheless, this work 
shows that a relativistic structure factor can be obtained with the help of the 
expressions derived in Refs. JU] an d by the use of simple mathematical proper- 
ties of the convolution integrals that describe the physical processes mentioned 
above. In both cases, extensive use has been made of the intracluster gas. 

To present these ideas we have divided the paper as follows. In section II, 
for reasons of clarity we summarize the ideas leading to the SZ spectrum for 
a non-relativistic electron gas emphasizing the concept of structure factor. In 
section III we discuss the relativistic case using an appropriate equation for 
the structure factor whose full derivation is still pending but leads to results 
obtained by other authors using much more complicated methods. Section IV 
is left for concluding remarks. 

2 General Background 

As it was clearly emphasized by the authors of this discovery in their early 
publications pQ [2] as well as by other authors, the distortion in the CMBR 
spectrum by the interaction of the photons with the electrons in the hot plasma 
filling the intergalactic space is due to the diffusion of the photons in the plasma 
which, when colliding with the isotropic distribution of a non-relativistic electron 
gas, generates a random walk. The kinetic equation used to describe this process 
was one first derived by Kompaneets back in 1956 [Z]-[H]- For the particular 
case of interest here, when the electron temperature Te ~ 10 8 A" is much larger 
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than that of the radiation, T = T^ad (2.726 K), the kinetic equation reads [3] 

dN_^ 2 d 2 N dN 
dy dv 2 dv 

Here N is the Bose factor, N — (e x — , x = fe, v is the frequency, h 

is Planck's constant, k Boltzmann's constant and y the "Compton parameter" 
given by, 

kTe f kTe , 
JJ= n cr T n e cdt= rT (2) 

m e being the electron mass, c the velocity of light, <7t the Thomson's scattering 
cross section and n e the electron number density. The integral in Eq. J2J) is 
usually referred to as the "optical depth, r", measuring essentially how far in 
the plasma can a photon travel before being captured (scattered) by an electron. 
Without stressing the important consequences of Eq. (JIJ readily available in 
many books and articles 0-0, we only want to state here, for the sake of 
future comparison, that since the observed value of N {v) is almost the same as 
its equilibrium value N eq (is), to a first approximation one can easily show that 



m = N{ V )-N eq {v) 
N ~ N eq {v) 



x 2 e x (e x + 1) Axe* 
(e x - l) 2 e x - 



(3) 



implying that, in the Ray leigh- Jeans region [x << 1), 



SN 

and in the Wien limit (x >> 1) : 

(5) 

If we now call I {v) the corresponding radiation flux for frequency defined as 

^) = ^ U v (T) (6) 

where U V {T) = 87r J 1/3 N eq (v) is the energy density for frequency v and temper- 
ature T , and noticing that 



ST _ f d(ln I{u))\ (SI 



T V dQjxT) J \ I 



(7) 



we have that the change in the background brightness temperature is given, in 
the two limits, by 

ST 

— = -2y 7 x«l (8) 



T 

and 



— =xy,x»l (9) 
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showing a decrease in the low frequency limit and an increase in the high fre- 
quency one. Finally, we remind the reader that the curves extracted from Eq. 
||3J| for reasonable values of the parameter "y" show a very good agreement with 
the observational data. 

As mentioned before, many authors, including Sunyaev and Zel'dovich them- 
selves, were very reluctant in accepting a diffusive mechanism as the underlying 
phenomena responsible for the spectrum distortion. The difficulties of using dif- 
fusion mechanisms to study the migration of photons in turbid media, specially 
thin media, have been thoroughly underlined in the literature 231" El- Sev- 
eral alternatives were discussed in a review article in 1980 [5] and a model was 
set forth by Sunyaev in the same year based on the idea that Compton 

scattering between photons and electrons induce a change in their frequency 
through the Doppler effect. Why this line of thought has not been pursued, 
or at least, not widely recognized, is hard to understand. Two years ago, one 
of us (ASV) 55 reconsidered the single scattering approach to study the SZ 
effect. The central idea in that paper is that in a dilute gas, the scattering law 
is given by what in statistical physics is known as the dynamic structure factor, 
denoted by G{k,v) where k = and A is the wavelength. In such a system, 

/ _—\2 

this turns out to be proportional to exp(— W V ) where W is the broadening 
of the spectral line centered at frequency v given by 

w= 2 2kTe 1/2v 

C TO 

If one computes the distorted spectrum through the convolution integral 

/>oo 

i»= / I {v)G{k,D-{l-ay)v)dV (11) 
Jo 

where the corresponding frequency shift ^7 has been introduced through Eqs. 
I|8I9|) (a = — 2 + x), one gets a good agreement with the observational data. This 
result is interesting from, at least, two facts. One, that such a simple procedure 
is in agreement with the diffusive picture. This poses interesting mathematical 
questions which will be analyzed elsewhere, specially since the exact solution to 
Eq. (1) is known (see Eq. (A-8) ref. [S]). The other one arises from the fact 
that this is what triggered the idea of reanalyzing the SZ effect using elementary 
arguments of statistical mechanics and constitutes the core of this paper. 

Indeed, we remind the reader that if an atom in an ideal gas moving say 
with speed u x in the x direction emits light of frequency v Q at some initial 
speed u x (0), the intensity of the spectral line / (v) is given by 

(12) 

Eq. (|12l) follows directly from the fact the velocity distribution function in 
an ideal gas is Maxwellian and that u x , u x (0) and v are related through the 
Doppler effect. For the case of a beam of photons of intensity I Q (v) incident 



loiy) 



exp 



2kT 



4 



on a hot electron gas regarded as an ideal gas in equilibrium at a temperature 
T e the full distorted spectrum may be computed from the convolution integral 
given by 



1 [°° 

1 M = Tun ^ / Io M exp 



f(v)v 



W(v) 



dv (13) 



which defines the joint probability of finding an electron scattering a photon 
with incoming frequency v and outgoing frequency v multiplied by the total 
number of incoming photons with frequency v. -j^W (v) is the normalizing 
factor of the Gaussian function for f(y) = 1, W (v) is the width of the spectral 
line at frequency v and its squared value follows from Ea. ljTU|) 

AhT 

W 2 (v) = ^rv 2 = Ayv 2 (14) 
vn e c 

where r is the optical depth whose presence in Eq. I|14fl will be discussed later. 
The function /(y) multiplying v in Eq. (|13[) is given by f(y) = 1 + ay where 
a = —2 in the Ray leigh- Jeans limit and a = xy in the Wien's limit, according 
to Eqs. (8-9) and the fact that ^ = for photons. Eq.JTSJ is the central 
object of this paper so it deserves a rather detailed examination. In the first 
place it is worth noticing that I {v), the incoming flux, is defined in Eq. ©. 
Secondly, it is important to examine the behavior of the full distorted spectrum 
in both the short and high frequency limits. In the low frequency limit, the 
Rayleigh- Jeans limit I {v) = 2k J/ , so that performing the integration with 
a = —2 and noticing that y is a very small number, one arrives at the result 

1151) is in complete agreement with the value obtained using Eq. (El ; t nc photon 
diffusion equation. In the high frequency limit where a = xy and / (^) = 
^%-e~ x , a slightly more tedious sequence of integrations leads also to a result 
at grips with the diffusion equation, namely 

51 

-j = x 2 y, x » 1 (16) 

Why both asymptotic results, the ones obtained with the diffusion equation and 
those obtained from Eq. i|13f) agree so well, still puzzles us. At this moment 
we will simply think of them as a mathematical coincidence. Nevertheless, it 
should be stressed that in his 1980 paper, Sunyaev [IJJj reached rather similar 
conclusions although with a much more sophisticated method, and less numer- 
ical accuracy for the full distortion curves. The distorted spectrum for the 
CMBR radiation may be easily obtained by numerical integration of Eq. I|13fl 
once the optical depth is fixed, y is determined through Eq.J5J) and a = —2 + x. 
The intergalactic gas cloud in clusters of galaxies has an optical depth r ~ 10~ 2 

Q- 
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From the results thus obtained in the non-relativistic case, one appreciates 
the rather encouraging agreement between the observational data and the the- 
oretical results obtained with the three methods, the diffusion equation, the 
structure factor or scattering law approach, and the Doppler effect. This, in our 
opinion is rather rewarding and some efforts are in progress to prove the math- 
ematical equivalence of the three approaches. From the physical point of view, 
and for reasons already given by many authors, we believe that the scattering 
Doppler effect picture does correspond more with reality, specially for reasons 
that will become clear in the last section. 



3 The relativist ic case 

Following the discussion of the previous section, Ea. H13l) may be written as: 

/>oo 

I(y)= / I {v)G{v,v)dv (17) 



\ W(u) ) 



and I n (v) is defined 



where the function G(y^ v) — ^w(u) ex P 

in Eq.(EJ). It is now clear that in Wien's limit, when jfe >> 1 , Eq. (|17J) becomes 
the Laplace transform of I a (y) with a parameter s = -g^ so that calling I p w 
the distorted spectrum in that limit, one gets that: 

oh f°° 

I pW (v) = — D 3 G(p, v) e~ sD dD (18) 
c Jo 

Comparing Eq. (|18|) with the result obtained by integrating the Kompaneets 
equation for the change of intensity in Wien's limit, namely [21-0 

oh 3 T 3 p x 

where F(x) = xcoth(f), a straight forward inversion of the Laplace transform 
yields that 



G(p, v) = S(p -v)- 4y^S'{D -u) + y%&" [y - u) (20) 

In this equation, the primes denote derivatives with respect to v . It may be 
thus regarded as a representation of G{v, v) in terms of the Delta function an 
its derivatives. To illustrate this point we can show that the RJ and Wien's 
limits of integral i|ll[) arise by substituting G(v, v) by a shifted delta function 
(see appendix). 

In the relativistic case, the form for G(D, v) may also be written down taking 

2 

the relativistic form for the particle's energy E = T "° c 2 and performing a 

V 



G 



2 jr rp 

power series expansion in powers of \. Defining z = ° 2 , this leads to an 
integral for I(v) which reads as follows: 



= TWi \ / Io i v ) exp 



I _ l*-f(y)vY 



2ry 



dv (21) 



Nevertheless, its integration even in the two limits Wien's and Rayleigh- 
Jeans has been unsuccesful. However, by direct numerical integration one 
obtains curves that agree qualitatively well with those obtained by other meth- 
ods. Thus, in order to verify the generosity of the delta function representation 
we proceeded in an indirect way. It has been shown in the literature that a good 
representation of the relativistic SZE up to second order in y reads as |1(J) : 

AI = y 2g£x 4 I -^[-4 + F (x) + if (-10 + %F{x) m 
-fF\x) + 1F 3 (!) - f H*(x) + lF(x)H 2 (x))} 

Here, H(x) — x [sinh(a?/2)] . Taking the inverse Laplace transform in 
Wien's limit of Eq. 

one obtains that 



G R (P, v) = Sip -v)- (Ay + i^) £8'(p - u)+ 

(v + ¥) - - -u) + ^ IV) (9 - v) 



(23) 



There is of course the remaining problem, namely to obtain Eq. I|22|l 
from the convolution integral (|21|l . This has so far defied our skills, but is purely 
a mathematical problem. The main point we want to underline here is the pos- 
sibility of using expressions such as Eqs. ( 120(1 and l(23|l . in principle obtainable 
from convolution integrals with structure factors which provide a useful ana- 
lytical tool that avoids resorting too complicated and laborious methods. In 
fact, Fig. 1 shows CMBR distortion comparison curves obtained using Eqs. ( 
I19|) and (|17|l with (|20(l . for several non-relativistic clusters. Fig. 2 shows the 
corresponding comparison relativistic curves at higher, realistic, electron tem- 
peratures using Eqs. ( 122(1 and ((17|l with 1(23(1 . It is clear that, for any practical 
purpose, the curves are identical. Indeed, one is not able to identify two different 
plots in each figure. 



4 Discussion of the results 

It is known that, for the case of a beam of photons of intensity I (y) incident 
on a hot electron gas regarded as an ideal gas in equilibrium at a temperature 
T e , the full non-relativistic SZE distorted spectrum may be computed from the 
convolution integral given in Eq. ((13(1 . which defines the joint probability of 
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Figure 1: CMBR SZE distortion as computed by Eqs. (5,10,12), with r = 1CT 2 , 
T = AKeV (lower curve), T — 6KeV (middle curve) and T — 8KeV (upper 
curve). No difference exists if compared the corresponding plot of Eq. (19). SI 
is measured in erg s _1 cm~ 2 ster^ 1 . The figure is scaled by a factor of 10 18 . 

I-Io 




Figure 2: Rclativistic CMBR SZE distortion as computed by Eqs. (5,11,13), 
with t = 1CT 2 , T = l2.5KeV (lower curve), T = IbKeV (middle curve) and 
T = 18KeV (upper curve). No difference exists if compared the corresponding 



plot of Eq. (22). SI is measured in erg s 
scaled by a factor of 10 18 . 



cm 2 ster . The figure is again 
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finding an electron scattering a photon with incoming frequency v and outgoing 
frequency v multiplied by the total number of incoming photons with frequency 
v. 1 — 2y is a function that corrects the outgoing photon frequency due to the 
thermal effect (u) is the normalizing factor of the Gaussian, W (u) is 

the width of the spectral line at frequency v and its squared value is defined in 
Eq. dJl. 

Thus, it is clear that an accurate convolution integral between the undis- 
torted Planckian and a regular function exists and can be obtained from strictly 
physical arguments in the non-relativistic SZE. It is interesting, however, that 
the corresponding relativistic expression does not seem to be a simple physical 
generalization of Eq. I|13|) . Yet, Eqs. (|17I23I) give the correct mathematical de- 
scription of the relativistic distortion by means of Dirac delta functions and its 
derivatives. Two questions can now be posed. One, of a rather mathematical 
fashion, states under what conditions a structure factor such as that appearing 
in Eq. I|17[l can be written in a representation involving Dirac's delta functions 
and its derivatives. This subject seems to be related to the mathematical the- 
ory of distributions. The second one, of more physical type, concerns with the 
possible description of the thermal SZE as a light scattering problem in which 
CMBR photons interact with an electron gas in a thermodynamical limit that 
allows simple expressions for the structure factors. If this were the case, the 
thermal relativistic SZE physics would not need of Montecarlo simulations or 
semi-analytic methods in its convolution integral description. 

Emphasis should be made on the fact that the main objective pursued in 
this work is to enhance the physical aspects of the SZE by avoiding complicated 
numerical techniques in many cases. This is achieved, as shown, by resorting 
the concept of structure factors, or scattering laws, widely used in statistical 
physics. 

APPENDIX 

The starting point here is Ea. l(T7|) . In the RJ limit, I oRJ (v) = ^£-D 2 . 
Introducing G(D 1 v) = 6(9 — v(l — ay)) one obtains 

IrM = -^r / v 2 5(v - v(\ - ay))dv = - ayf 

o 



keeping terms linear in y we find that 

Irj(v) - IoRJ (v) 



loRJ 



-2ay 



This result is consistent with Eq. (8) setting a = 1. 

Now, in the Wien limit, I Q w (v) — fr^e - ^ , in this case the introduction 
of 8(i> — v(l — ay)) leads to 

I w (i/) = ^ ( v 3 e~&S(D - v(l - ay))dD =^v 3 (l- ay) 3 e~&e£r ay 
c Jo c 
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Expanding e a v arL d keeping terms up to first order in y we obtain that 

Iw{v) - hw (y) a 

r 7 \ = -6ay + axy 

low \y) 

Finally, setting a = 1 and considering x >> 6 in the Wien limit, we obtain the 
desired result, Eq. 

This work has been supported by CONACyT (Mexico), project 41081-F. 
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